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Testing the local-void alternative to dark energy using 
galaxy pairs 
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ABSTRACT 

The possibility that we live in a special place in the universe, close to the center of a 
large, radially inhomogeneous void, has attracted attention recently as an alternative 
to dark energy or modified gravity to explain the accelerating universe. We show that 
the distribution of orientations of galaxy pairs can be used to test the Copernican 
principle that we are not in a central or special region of Universe. The popular void 
models can not fit both the latest type la supernova, cosmic microwave background 
data and the distribution of orientations of galaxy pairs simultaneously. Our results 
rule out the void models at the 4a confidence level as the origin of cosmic acceleration 
and favor the Copernican principle. 

Key words: cosmology: theory - dark energy 



> 

ON 

o> 
m 

o 



* 



1 INTRODUCTION 

The standard model of cosmology based on the cosmological 
principle (homogeneity, isotropy, validity of General Rela- 
tivity) which contains about 23% dark matter, 4% ordinary 
matter and 73% dark energy driving the acceleration of a 
flat universe has been established. Many astronomical ob- 
servations support this standard picture, including type la 
supernovae (SNe la) (Riess et al. 1998; Perlmutter et al. 
1999), cosmic microwave background (CMB) (Komatsu et 
al. 2011; Sherwin et al. 2011), baryon acoustic oscillations 
(BAO) (Eisenstein et al. 2005) and gamma-ray bursts (Dai 
et al. 2004; Wang et al. 2007; 2011). 

In the meanwhile, inhomogeneous Lemaitre-Tolman- 
Bondi (LTB) (Lemaitre 1933; Tolman 1934; Bondi 1947) 
universe could also induce an apparent dimming of the light 
of distant supernovae. The idea is to drop the dark en- 
ergy and the Copernican principle, and instead suppose that 
we are near the center of a large, nonlinear ly underdense, 
nearly spherical void surrounded by a flat, matter domi- 
nated Einstein-de Sitter (EdS) spacetime. Because the ob- 
server must be at the center of void, so the LTB models 
violate the Copernican principle. Because of the observed 
isotropy of the CMB, the observer must be located very close 
to the center of the void (Alnes &Amarzguioui 2006). It was 
demonstrated LTB models can fit the SNe la data, as well as 
the BAO data and the CMB data (Garcia-Bellido & Haug- 
boelle 2008). Some tests have actually been proposed: the 
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Goodman-Caldwell-Stebbins test, which looks at the CMB 
inside our past lightcone (Goodman 1995; Caldwell & Steb- 
bins 2008), the curvature test, which is based on the tight 
relation between curvature and expansion history in a Fried- 
mann spacetime (Clarkson et al. 2008), and the radial and 
transverse BAO scale (Zibin et al. 2008; Garcia-Bellido & 
Haugboelle 2009). However, based on these tests, void mod- 
els have not yet been ruled out (Clifton et al. 2008; Uzan et 
al. 2008; Biswas et al. 2010; Wang & Zhang 2012; Nadathur 
& Sarkar 2011). Zhang & Stebbins (2011) have excluded 
the Hubble bubble model as the possibility of cosmic accel- 
eration using the the Compton-y distortion. Zibin & Moss 
(2011) also concluded that a very large class of void models 
was ruled out using this method. Here we propose a powerful 
tool, the orientations of galaxy pairs to test the Copernican 
principle. 

The Alcock-Paczynski (AP) test is a purely geometric 
test of the expansion of the Universe (Alcock & Paczyski 
1979). Marinoni & Buzzi (2010) implemented the AP test 
with the distribution of orientations of galaxy pairs in orbit 
around each other in binary systems. The principle of this 
method is that the orientations is thought to be completely 
random, with all orientations being equally likely if mea- 
sured assuming a cosmology that matches the true under- 
lying cosmology of the Universe in a Friedmann-Lemaitre- 
Robertson- Walker (FLRW) universe after the effect of pe- 
culiar motion is excluded. 

In this paper, we implement the Alcock-Paczynski test 
with pairs of galaxies to test the Copernican principle. The 
void models cannot both fit SNe la plus CMB data and ori- 
entations of galaxy pairs. Our results exclude the possibility 
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of the void models as the source of cosmic accelerating ex- 
pansion and favor the Copcrnican principle. 



2 THE VOID MODEL 

We model the void as an isotropic, radially inhomogcncous 
universe described by the LTB metric, 
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where a prime denotes the partial derivative with respect to 
the coordinate distance r, and the curvature k(r) is a free 
function representing the local curvature. The transverse ex- 
pansion rate is defined as H± = A(r, t)/A(r, i) and the radial 
expansion rate is defined as Hu = A'(r,t)/A'(r,t), where an 
overdot denotes the partial derivative with respect to t. 

The Friedmann equation in LTB metric is H\ — 
F(r)/A 3 (r,t) + c 2 k(r)/A 2 (r,t), where F(r) > is a free 
function which determines the local energy density. The di- 
mensionless density parameters can be determined as Qm (r) 
and n K (r) by F(r) = Hl{r)Sl M {r)Al(r) and c 2 k{r) = 
Hf ) {r)Q,K{r)A 2 ) {r), where Ho(r) and Ao(r) are the values of 
H±(r,t) and A(r,t) respectively at the present time t = to- 
So we can rewrite Friedmann equation in LTB metric as 
Hl(r,t) = Hi(r)[Q M (r)(A /A) 3 + Q K {r)(A /A) 2 ]. This 
equation can be integrated from the time of the Big Bang, 
£b = £b(V)> to yield the age of the universe at any given 
(r,t), 
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The function, Ao(r), corresponds to a gauge mode and we 
choose to set Ao(r) = r. As stressed by Silk (1977) and Zibin 
(2008), it is crucial to consider only voids with vanishing 
decaying mode, so we set £b(t") = everywhere. Although 
Biswas et al (2010) have shown that the void models were 
in better agreement with observations if the void has been 
generated sometime in the early universe. The null radial 
geodesies described by 
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where Hu(z) — Hu (r(z),t(z)). The angular diameter dis- 
tance and luminosity distance are given by 

d A (z) = A (r(z), t(z)) , d L (z) = (1 + zf A (r(z), t(z)) . (4) 

We will adopt the two parameterizations of the void 
profile Qjuir)- The first one is the constrained GBH model 
(Garcia-Bellido & Haugboelle 2008) 



^m{t) — ^iM,out+(^M,in — ^M,out) ■ 



- tanh(r - r /2Ar) 
1 + tanh(r /2Ar) 



(5) 



where the parameters r\) and Ar characterize size and steep- 
ness of the density profile respectively. We wish to look only 
at voids that are asymptotically EdS, so we set Qm, out = 1- 
We also set Ar = 0.35ro, because this value can well fit 
the SNe la data (Garcia-Bellido & Haugboelle 2008; Marra 
& Paakkonen 2010). This density shape of GBH model can 
also explain other observations, such as CMB and BAO. The 
second one is a simple Gaussian form, 
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Figure 1. The la, 2a and 3a contours in the parameter space 
&M,in — fo f° r the constrained GBH model. The dash-dot contours 
represent constraint from SNe Ia+CMB, the dashed contours 
from AAP, and the solid contours from SNe Ia+CMB+AAP. The 
best fit parameters arc tq = 3.0 Gpc and S!m j n = 0.10 for SNe 
Ia+CMB. But in order to fit the AAP, much more larger and 
undcrdense void arc needed. This void model can not fit both the 
SNe la plus CMB data and the orientations of galaxy pairs. 



n M (r) = M ,out + (n M ,in — n M ,out) exp(-r /r ), 



(6) 



where Qmm and flM,out are the matter density parameters 
at the observer's position and in the FLRW background out- 
side the void, and ro characterizes the size of the void. It has 
been shown that this void profile can fit the observations of 
SNe la, CMB and BAO (Nadathur & Sarkar 2011). 



3 CONSTRAINT FROM SNE IA, CMB AND 
PAIRS OF GALAXIES 

We use the recent Union2 SNe Compilation (Amanullah et 
al. 2010), which consists 557 SNe la in the redshift range z — 
0.015 — 1.4. With dh in units of megaparsecs, the predicted 
distance modulus is jj,(z) — 5 log d^ (z) + 25. The likelihood 
analysis is based on the \ 2 function: 
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The parameter jj, is an unknown offset. We marginalize the 
likelihood exp(— Xs 2 No/2) over /i, leading to a new marginal- 
ized \ 2 function: 

S 2 

XSNe = $2 - -g- , (8) 

where S„ - ^2 t [n(zi) - ^ob s {zi)] n / a 2 . 

We also use positions and amplitudes of peaks and 
troughs in the CMB spectrum to test the LTB models. The 
location of peaks and troughs can be calculated as Hu et al. 



d A (z*)(l + z") 



where 



(2001): Im — (m — 4> m ) Ia , where Ia — n 
dA(z*) is the angular diameter distance with the sound hori- 
zon of r* at the recombination redshift of z* . We use the 
method of Marra & Paakkonen (2010) to calculated these 
values. We consider the position of the first, second, third 
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Figure 2. Same as Fig. 2, but for the gaussian LTB model. The 
best fit parameters are ro = 3.6 Gpc and S7M,in = 0.12 from SNe 
Ia+CMB. The AAP favors a much more larger and underdense 
void. 



peak and of the first trough. We compute the correspond- 
ing phases <j>\, cf>i.5, 02 and (f>3 using the accurate analytical 
fits of Doran & Lilley (2002) . The relative heights of second 
and third peak relative to the first one, Hi and Hz are also 
considered, for which we can use the fits of Hu et al. (2001). 
So the Xcmb i s (Marra & Paakkonen 2010) 
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where the W7 represents the best-fit WMAP7 spectrum 
(Jarosik et al. 2011). 

Pairs of galaxies should be distributed with random ori- 
entations if the fundamental assumptions of homogeneity 
and isotropy are correct. But two factors affect this simple 
cosmology test. First, peculiar velocities displace the posi- 
tion of a galaxy along the line of sight from its true position. 
Marinoni and Buzzi modelled the peculiar velocity distortion 
as a Doppler shift where the observed line of sight separa- 
tion is related to the actual separation. Second, an observer 
needs to assume a cosmological model to convert observed 
angles and redshifts into comoving distances. The uniform 
distribution of orientations is distorted if a wrong underlying 
cosmology of the Universe is assumed. 

In a non-flat ACDM universe, the tilting angle t sub- 
tended between galaxy pairs and the line of sight, can be 
written as 

sin 2 t = {1 + [C k ( XA )cot9 - gfc _ ( f ^f*^ ] 2 }- 1 , (10) 

for details, see Marinoni et al. (2012). It is nontrivial to 
calculate the tilting angle and average anisotropy of pairs 
in LTB models. The measured galaxy (matter) clustering 
and its evolution agree with the standard ACDM cosmol- 
ogy to a factor of about 2 uncertainty up to z~l (Tegmark 
et al. 2004; Coil et al. 2006; Fu et al. 2008; Schrabback et 
al. 2010; Guzzo et al. 2008). A minimalist approach is to 
simply use the ACDM value since any viable LTB models 
must be consistent with these data. So we use the observed 
average anisotropy of pairs from Marinoni & Buzzi (2010) 
derived in ACDM cosmology. Zhang & Stebbins (2011) also 



approximated the matter power spectrum by its form in a 
standard ACDM cosmology. The observed tilting angle is 
shifted to apparent angle r because of the geometric distor- 
tions induced by the peculiar velocities of the pair's mem- 
bers. The probability distribution function of the apparent 
angle r which is is given by Marinoni & Buzzi (2010) 

t/ si 1 (l + o- 2 )(l + tan 2 r) . ., .„„. 

*Mdr = - / \ ^- — -Utanr|dT, (11) 
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and the parameter a depends on the cosmological expansion 
history as 



o(z) — a 



H (r){l + z) 



(12) 



The normalization parameter a is given by a = 
Hq 1 ((dv 2 /dr 2 \) . Because in the LTB metric, the trans- 
verse and radial expansion rates are different, so the correct 
value must be used in our calculations. When we use the 
galaxy pairs for AP test, the velocity perturbation a used in 
equation (12) is related to the peculiar motions of the pair 
members along the line of sight. So we use radial expansion 
rate Hu to calculate the velocity perturbation. Because the 
AP test is similar to BAO, we can see this formula is also 
similar to the redshift interval 5z corresponding to the acous- 
tic scale in the radial direction (Garcia-Bellido & Haugboelle 
2008; Biswas et al. 2010; Marra & Paakkonen 2010). In the 
homogeneous ACDM model, Marinoni & Buzzi (2010) used 
the normal expansion rate H(z). Marinoni & Buzzi (2010) 
derived the distribution ^(r) as the average anisotropy of 
pair (AAP), which is given by 



fi a = / sin r^(r)dT = 



(1 + a 2 ) arctan(o-) 



(13) 



At z » 0, Marinoni & Buzzi (2010) obtained a = 5.79+q jj 2 , 
using binaries in the seventh data release of the Sloan Digi- 
tal Sky Survey (SDSS) (Abazajian et al. 2009). The normal- 
ization factor a is assumed to be constant for all redshifts 
and for different galaxy selections (Marinoni & Buzzi 2010). 
Although Jennings et al. (2012) found that the value of a 
could have a small variation with cosmology and redshift, 
Marinoni & Buzzi established that the changes of best fit 
value cannot exceed the la confidence level if the variation 
of a is less than 10%. So this assumption could be reason- 
able. Belloso et al. (2012) also found that observations of 
close-pairs of galaxies do show promise for AP cosmologi- 
cal measurements, especially for low mass, isolated galaxies. 
The high-redshift (up to z ~ 1.45) AAP are obtained using 
the third data release of the DEEP2 survey (Davis et al. 
2007). The value of Xaap is 
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(14) 



We adopt the value of Hcr,obs and a b s from Fig. 2 of Marinoni 
& Buzzi (2010), which are shown as points in the Fig [3] In 
order to verify the hypothesis that the normalization factor 
a is constant for all redshifts, the distance between the ob- 
served value of recession velocity difference square (dV 2 (z)) 
and the prediction of equation (S20) in Marinoni & Buzzi 
(2010) is minimal (see Marinoni & Buzzi (2010) for more 
details). So this x 2 value is 
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We use the value of (dVo(z)} and a dv i from Fig.(5S) of 
Marinoni & Buzzi (2010). The total xaap is 



2 _ l% i '2 

Xaap — Xaap + Xdv 2 • 



(16) 



In Fig. [T] we show the l<r, 2a and 3cr contours in the 
fW.in — ro plane for the constrained GBH model. In the 
calculation, the priors from WMAP7, such as the age of 
Universe to = 13.79 Gyr and spectral index n s — 0.96 are 
used (Komatsu, et. al. 2011). We also marginalize the Hub- 
ble constant Ho in the range 50 ^ Ho ^ 80 km s _1 Mpc~ . 
The constraint from SNe la and CMB is shown as dash-dot 
contours, and dashed contours for AAP. The allowed range 
of r is 1.80 Gpc < r < 4.10 Gpc at 3<r level from SNe 
Ia+CMB. But the allowed range of ro is ro > 4.42 Gpc at 
3er level from AAP. These two contours do not overlap. So 
the constrained GBH model can not explain the observations 
of SNe Ia+CMB and AAP. The solid contours are derived 
from SNe Ia+CMB+AAP with xLn = 641.40. While for the 
ACDM model, the minimum x 2 is 620.37. The constrained 
GBH model is excluded at the 4<r confidence level compared 
to ACDM. In Fig[5] we show the la, 2a and 3a contours 
in the JlM.in — ro plane for the gaussian LTB model. The 
solid contours are derived from SNe Ia+CMB+AAP with 



Xn 



= 646.50. This model is also excluded at the 4a confi- 



dence level compared to ACDM. From the Xmin °f t fle two 
void models, we conclude that the precise form of the density 
profile may not be essential. Because the void models depend 
crucially on the void depth fo = (fiM,in — £lM,out)/&M,out 
and the void size ro . So our conclusion is almost independent 
of void model. 

In Fig[3] we show the theoretical redshift scaling of the 
AAP in these two LTB models. In the up panel, we use the 
best fit parameters from SNe Ia+CMB for the constrained 
GBH model, ro = 3.0 Gpc and Qm, in = 0.10. Obviously, 
the predicted values of AAP deviate from the observational 
values at high redshift. The x 2 value is 37.35 for these nine 
data points. In the bottom panel, ro = 3.6 Gpc and Qm.iu = 
0.12 are used for the Gaussian LTB model. The x 2 value is 
41.96 for these nine data points. 

We must note that the local Hubble constant iifi oc is 
also a big obstacle to the void models. Because the measure- 
ment of the Hubble constant is carried out mostly within a 
distance of roughly r\ oc ~ 200 Mpc (Riess et al. 2011; Freed- 
man et al. 2012), we obtain the Hi oc (Marra & Paakkonen 
2010) 



H h 



H (r)4-Kr 2 dr/(4TT/3r? o , 



(17) 



In order to fit both the SNe la and CMB, the value of H\ oc is 
64±3.2 km s~ 1 Mpc~ 1 in the constrained GBH model or 63+ 
3.5 km s~ 1 Mpc~ 1 in the Gaussian LTB model. Riess et al. 
(2011) determined the Hubble constant with 3% uncertainty 
as 73.8+2.4 km s" 1 Mpc -1 . Freedman et al. (2012) measured 
the Hubble constant as 74.3 ± 2.1 km s" 1 Mpc -1 . 
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Figure 3. Points represent the observed average anisotropy of 
pairs. Solid lines represent the theoretical redshift scaling of the 
AAP as predicted by Eq. 1 1131) in different LTB models with best 
fit parameters from SNe Ia+CMB, up panel for constrained GBH 
and bottom panel for gaussian LTB. The dashed line shows the 
best fit ACDM model with U M = 0.25 and Q A =0.65 



4 DISCUSSIONS 

Previous investigations show that void models can fit a va- 
riety of cosmological observations without containing dark 
energy because the lack of homogeneity gives a great de- 
gree of flexibility. For example, since the last scattering sur- 
face is far away from regions where SNe la are observed, 
the property of inhomogeneity allows a model to be con- 
structed which provides different physical densities in the 
regions from which these two sets of observational data are 
drawn. So, the best way to constrain inhomogeneous models 
is using several sets of data that measure a range of observ- 
ables at comparable redshifts. In this paper, we confront two 
general classes of void models with observations of SNe la, 
CMB and orientations of galaxy pairs. The redshifts of SNe 
la and orientations of galaxy pairs are almost in the same 
range. We find that the these two void profiles can not fit 
both SNe la plus CMB data and orientations of galaxy pairs 
simultaneously. We also show that the two void models can 
fit both SNe la and CMB data, but at the expense of a Hub- 
ble constant so low that they can also be ruled out. So our 
results favor the Copernician principle. We must also note 
that our results are obtained under some assumptions, such 
as the chosen priors and void profile, which is also discussed 
in Biswas et al. (2010). So the void models are ruled out at 
the 4cr confidence level given the explored models and pri- 
ors. But observations challenge the void models (Biswas et 
al. 2010; Zibin & Moss 2011). Future galaxy surveys such as 
BigBOSS (Schlegel et al. 2011) will provide improved preci- 
sion of AAP function, placing much more strong constraints 
on inhomogeneity. 
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